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Abstract 

Recently the ATLAS collaboration reported a 3 cj excess in the search for the events contain¬ 
ing a dilepton pair from a Z boson and large missing transverse energy. Although the excess 
is not sufficiently significant yet, it is quite tempting to explain this excess by a well-motivated 
model beyond the standard model. In this paper we study a possibility of the minimal su¬ 
persymmetric standard model (MSSM) for this excess. Especially, we focus on the MSSM 
spectrum where the sfermions are heavier than the gauginos and Higgsinos. We show that the 
excess can be explained by the reasonable MSSM mass spectrum. 



1 Introduction 


Recently the ATLAS collaboration has reported a 3cr excess in the search for the events with a 
dilepton pair on Z boson mass peak and large missing transverse energy (MET) [1]. The signal 
events should contain a same-flavor opposite-sign dilepton pair with its invariant mass in the Z 
boson mass range, large MET (E™®® > 225 GeV), and large scalar sum (iLx > 600 GeV) of the 
transverse momenta of all signal jets and the two leading leptons. The observed events in the 
dielectron and dimuon channels are 16 and 13 respectively, whereas the expected numbers of the 
standard model (SM) background events are 4.2 ± 1.6 and 6.4 ± 2.2. The two channels combine to 
give 29 observed events compared to 10.6 ± 3.2 expected from the SM, amounting to a 3 ct excess. 

On the other hand, the GMS collaboration has an analogous search for events with large MET 
and a dilepton pair on Z [2], where the signal events are classihed into MET bins = 100-200, 

200-300 and > 300 GeV) and the number of jets njet > 2, 3. In this GMS counterpart search, no 
signihcant deviation from the SM expectation is found. 

Although the present ATLAS Z excess is not so statistically signihcant yet and might even 
conhict with the aforementioned counterpart search by the GMS collaboration, it is quite tempting 
to investigate whether a well-motivated model beyond the SM can explain the signal. The large 
MET is a typical signature of the supersymmetric (SUSY) SMs. Several studies are devoted to the 
ATLAS Z excess in the context of the SUSY SMs. 

In Ref. [3], general requirements in the SUSY SMs to explain the excess are studied, and it 
is concluded that there needs to be a particle lighter than about 1.2 TeV, with production cross 
section the order of colored particles, as well as producing 0(1) Z bosons in its decay chain. In the 
specihc case of the minimal SUSY standard model (MSSM) with a neutralino as the lightest SUSY 
particle (LSP), it is found that Z bosons are generically not produced enough in the SUSY cascade 
decay chain. One can try to compensate this by increasing the production cross section with the 
light gluinos or squarks, but the constraints from the jets and MET searches [4] get more severe. 
Possible solutions to overcome this difficulty are discussed, including use of a compressed spectrum 
or a spectrum with the light gravitino LSP, in which the neutralino next-to-LSP (NLSP) decays 
into the gravitino and Z boson. Another type of spectrum is studied in Ref. [5], where the hrst and 
second generation squarks decay into a Bino followed by the Bino decaying into Higgsinos and W, 
Z boson, or the Higgs boson h. 

The spectrum with the light gravitino LSP is realized in the so-called general gauge mediation 
(GGM) [6, 7], where the branching fraction of neutralino NLSP into Z can be close to one. However, 
constraints from other SUSY searches such as jets-|-MET [4], stop search [8] or multi-lepton [9, 10] 
as well as GMS on-Z dilepton [2] are severe and such an explanation is not viable [11]. 

On the other hand, compressed spectra are utilized in other attempts to explain the excess [12- 
14]. When the mass difference between the neutralino LSP and the neutralino NLSP is less than the 
Higgs mass (125 GeV), the two-body decay of the NLSP into the LSP plus a Z boson can be efficient. 
An important requirement here is to ensure that the parent particle (gluino or squark) decays mainly 
into the NLSP so that the NLSP can produce Z in the next step of the decay chain. In the case of the 
GGM with a light gravitino, this requirement is satished because the gravitino LSP is very weakly 
coupled. In the heavy LSP scenarios, the LSP is taken to be a Bino-like neutralino in the MSSM 
with light sbottom [13] or singlino-like in the next-to-MSSM (NMSSM) [12, 14]. Gombining other 
constraints [2, 4, 15], the NMSSM scenario can reduce the signihcance of the excess only in a small 
region in the parameter space [14]. The light sbottom scenario [13] can also reduce the signihcance. 
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Figure 1: Schematic picture of the present MSSM spectrum. 


but sbottom produces bottom quarks when it decays into Higgsino-like neutralinos, and hence is 
severely constrained by 6-jets searches [19]. There is also a non-SUSY study [16] in the composite 
Higgs [17] / Randall-Sundrum [18] framework, but this also has bottom-rich signatures. 

In this paper, we revisit the possibility of explaining the excess in the MSSM, in particular in 
the well-motivated split SUSY-like spectrum. The MSSM is one of the most attractive candidates of 
models beyond the standard model. Especially the recent discovery of the Higgs boson h with a mass 
of around 125 GeV [20] seems to suggest the framework of the split SUSY [21, 22], where the SUSY 
fermions are around TeV scale and the SUSY scalars are heavier than TeV scale. This framework can 
overcome weak points of the weak-scale MSSM, such as SUSY flavor/CP and cosmological problems. 
Most importantly, the split SUSY is quite compatible with the observed Higgs mass [23, 24]. In the 
light of the Higgs discovery, this framework is intensively studied [25-28]. Now the MSSM with split 
SUSY-like spectrum is one of the most convincing and viable models. Therefore it is very interesting 
to study this model in light of the ATLAS Z excess. 

To explain the ATLAS Z excess, the MSSM mass spectrum will at least satisfy the condition 
that the SUSY cascade decay chain has large branching fraction to Z. Even if the Z-rich decay 
chains are realized at the LHC, the non-leptonic decays of the Z bosons induce signals of multi-jets 
plus MET with zero lepton, which are severely constrained. Moreover, multiple leptonic decays of 
the Z bosons may result in multi-lepton events. Therefore our goal is to hnd the mass spectra which 
satisfy the conditions: 

• SUSY cascade decay chain is Z-rich. 

• Less constrained by searches other than the dilepton channel, such as multi-jets-fMET. 

In this paper, we point out that a simple mass spectrum (Fig. 1) can well account for the ATLAS 
Z excess. If the mass spectrum is compressed enough, ~ 100 GeV and > rrig — 2mt, 

the above conditions can be simultaneously satished. As we will see in Section 2, when the stop is 
heavy and the spectrum is sufficiently compressed, the gluino radiative decay into a Higgsino-like 
neutralino and a gluon becomes the dominant gluino decay channel. The Higgsino-like neutralino 
can then decay into Z boson with a branching fraction close to 1. Together, this gives rise to an 
efficient production of Z boson. In Section 3, we reduce essential features of the MSSM spectrum to a 
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simplified model, and study LHC signals of the model and its constraints. Summary and discussions 
are given in Section 4. 

2 SUSY Spectrum 

2.1 Mass spectrum 

We take the split SUSY-like spectrum, where the gauginos and Higgsinos are light whereas scalar 
superparticles are heavy, and we consider production of the relatively light gluino, which decays into 
neutralinos and charginos. The LSP is a Bino-like neutralino and there are nearly degenerate 
two Higgsino-like neutralinos X 2 3 a Higgsino-like chargino xf cis the NLSPs. For simplicity, we 
take the Wino heavier than the gluino. 

Because of the renormalization group effects, the right-handed stop is expected to be typically 
lighter than the other squarks. If the right-handed stop is lighter among the squarks, gluino branch¬ 
ing fraction shows an interesting feature: the dominant decay channel becomes the top-stop-loop- 
induced process into a Higgsino-like neutralino and a gluon for suitable mass splitting between gluino 
and the neutralinos [29-31] (see Fig. 2). This reduces the branching fraction of gluino three-body 
decay modes, and the gluino dominantly decays into a neutral Higgsino with a gluon. 



Figure 2: Diagrams of the gluino decay into a Higgsino and gluon. 


2.2 Decay of neutralinos and gluino 

We study neutralino and gluino decays in this Subsection to motivate the present split SUSY spec¬ 
trum and a simplihed model studied in the next Section. 

Decay of neutralinos 

First, we consider the neutralino decay. A neutralino may decay into a lighter neutralino emitting a 
Z or h boson, and into a chargino with a W boson if each channel is kinematically allowed. In the 
present mass spectrum, the lightest chargino xf is Higgsino-like, so it is approximately degenerate 
with Higgsino-like neutralinos X 2 3- This means that neutralinos cannot decay into the chargino with 
a W boson in the present spectrum. If the mass difference between Xi and ^ 3 is greater than mz 
and less than rrih, BF(x 2 3 Xi-^) — 1- Once the Higgs channel becomes kinematically available, 

the branching ratio of each Higgsino-like neutralino into Z or h varies substantially in the case of 
low tau/d depending on the phases of parameters such as /i-term, but X]i =2 3 ^^(Xi' x^Z)/2 ~ 
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Xi^)/2 — 0.5 in the limit of large mass difference between Xi and X 2 3 - This can 
be understood in the Nambu-Goldstone picture: the longitudinal component of the Z boson is the 
complex partner of the Higgs, whereas the transverse components are unimportant in the limit. 
Taking the average is justihed because gluino decays into the up-type Higgsino in our spectrum 
(see below), and it is approximately equally contained in the two mass eigenstates of Higgsino- 
like neutralinos. Therefore in order to produce the Z bosons in the gluino decay chain efficiently, 
mz < m^o ^ — m^o < rrih is required. In the following analysis, we assume m^o ^ — m^o ~ 100 GeV. 


Decay of gluino 

Next, let us move on to the gluino decay. As studied in Refs. [29-31], the partial decay rate of gluino 
into a gluon and an (up-type) Higgsino is relatively enhanced by a factor (log {mi/rrit))^ compared 
to other channels: (i) gluon and Bino, (ii) neutralino, quark and antiquark, and (iii) chargino, quark 
and antiquark. Therefore, the gluino efficiently produces Higgsino-like neutralinos X 2 X 3 the 
case of the light gluino and heavy (but relatively lighter among the squarks) stop. This radiative 
decay of the gluino has some advantages to explain the ATLAS Z excess. The gluino decay into a 
Higgsino-like neutralino, followed by the decay of the neutralino into a Z boson, efficiently produces 
the Z boson. Thanks to the log enhancement, this decay mode dominates over the other channels 
for wide parameter space and the expectation value of the number of Z bosons per gluino decay is 
enhanced. In this case, the constraints from other SUSY searches, such as multi-jets and leptons 
signals, can be relaxed. Another advantage of this radiative decay is that we can suppress the 
branching fractions into heavy-flavor jets, which are severely constrained by LHG searches even for 
the compressed mass spectrum. 

In the estimation of gluino branching fractions, we need resummation of the log {m^/rrit) factor, 
otherwise, the two body decay rate is overestimated [30] . For this resummation, we first evaluate the 
Wilson coefficients of the dimension five dipole operator Oj = dimension six dipole 

operator and several four-Fermi operators which include a gluino spinor e.g., 

at the sfermion mass scale. The operators relevant for the radiative 

decay of the gluino are Of, and At the sfermion mass scale, these Wilson coefficients 

are given by 


'-' 2,33 
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Here we neglect sfermion flavor violation and the Yukawa couplings other than yt. By using the 
renormalization group equations (RGBs) of the Wilson coefficients, we evolve down to the gluino 
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mass scale. The RGEs of interest are written as, 

(iln/x ^ IGtt^ ^ ’ 

^ (c,%\ + ^9syt \ (c^,h\ 

rfln/x j Ihvr^ \ AgsVt -lAgl + ^yj) j ' 


( 4 ) 

( 5 ) 


Here we take into account only the top Yukawa and strong couplings. The Wilson coefficient (7^“ is 
logarithmically enhanced through the operator mixing with 0^33 and Cf is not. Thus, the radiative 
decay rate of the gluino into the Higgsino are relatively enhanced with heavier sfermions. We 
evaluate the gluino decay rates with these Wilson coefficients at the gluino mass scale. For details, 
see Refs. [30, 31]. Following Ref. [32], we set the weak scale Yukawa and gauge coupling parameters. 

In the calculation of gluino branching fractions, the relevant parameters are gluino mass rrig, 
Bino and Wino mass parameters Mi and M 2 , squark masses, Higgs fi parameter, CP odd Higgs 
mass niA, and tan/1. For simplicity, we neglect mixing among squarks, and we take a universal mass 
for squarks except for the right-handed stop. The former is set three times larger than the latter, 
and rrig = 800 GeV is taken unless otherwise stated. To reduce the number of parameters, we fix 
M 2 = 3 TeV, and uia same as squark masses. To make mass difference between Xi and X 23 ) xf 
about 100 GeV, we £x |/i| = Mi -|- 100 GeV with Mi positive. We vary Mi and diagonalize the 
neutralino (and chargino) mass matrix to show the xtilsp axis in the Figure 3, where the branching 
fractions of gluino are plotted. 

In Fig. 3a, hve dominant branching fractions of gluino are shown with the right-handed stop mass 
= 10 TeV. The black line represents the branching fraction of the gluino to the Higgsino-like 
neutralinos and gluon. It diminishes in the right side of the Figure for the kinematical reason. As 
the mass splittings between the relevant neutralino or chargino and gluino increase (to the left of 
the Figure), the three-body decay channels emitting a quark-antiquark pair become non-negligible. 
The choice of the parameters in Fig. 3a is relatively inefficient for the gluon channel. The branching 
fraction of the gluino into gxts other choices of parameters are shown in Fig. 3b. The branching 
fraction becomes large in particular when we take the relative sign between /i and Mi negative. This 
is because the up-type Higgsino component in the LSP decreases by partial cancellation in the case 
of low tan/3 and sgn(/x/Mi) = —1. 

In Fig. 3c, we show the dependence of the gluino branching fraction of ^ gX 2 3 ^^e squark 
masses. Here we assume rrig^ 2 = rriq^^ = m^. Mi = 500 GeV, /i = —600 GeV and tan (3 = 2. 
We show the cases that rriq/m^^ = 1,1.5,2 and 3. If the rriq/mi^ is large enough, = 0(1) TeV 
can realize the favoured gluino decay, which may leave the possibility that the relatively “natural” 
SUSY can account for the ATLAS Z excess. As seen in Fig. 3c, heavier stop can increase the 
BF(^ —( 7 X 23 ), ‘iue to the large log enhancement. Note that, however, the gluino decay length gets 
larger as the stop mass increases: 


CTg ~ O(10)/im 


frrig — ruNLSP 

V 300 GeV 


100 TeV 
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( 6 ) 


The standard tracking system assumes the gluino decay occurs within 0(1) mm from the primary 
vertex. Therefore the stop mass should be less than 0(100) TeV to produce the standard MET 
and/or jets and/or leptons signals. If the decay length is larger than 0(1) mm, severer constraints 
will be imposed even in the case of the compressed mass spectrum [33]. 
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(a) Branching fraction of each decay mode (b) LSP mass dependence of BF(^ —)• 



Figure 3: Various dependencies of branching fractions of the gluino. (a): The branching frac¬ 
tions of the gluino decay as functions of the LSP mass utlsp = m^o. The solid black, long- 
dashed red, medium-dashed blue, short-dashed green, and dot-dashed pink lines show the branch¬ 
ing fractions to ( 7 X 2,35 5 9Xii ttxi, and tbxf, respectively. We take rrig = 800 GeV, and 

(sgn(/i), tau/S,m^^/TeV) = (-[-,2,10). (b): The branching fraction into ( 7 X 23 various choices of 
parameters. The solid black, long-dashed red, medium-dashed blue, and short-dashed green lines 
correspond to (sgn(/r), tan/3, m^^/TeV) = (—, 2,10), (—, 50,10), (—, 2,100), and (-|-,2,100), respec¬ 
tively. (c): The branching fraction into 77 X 23 ^ function of the stop mass The solid black, 

long-dashed red, short-dashed blue, and dotted green lines correspond to squarks masses 1, 1.5, 2, 
and 3 times larger than the lighter stop mass. We set Mi = 500 GeV, /i = —600 GeV, and tan/3 = 2. 
(d): The branching fraction of ^ 77 X 23 fh® ^3 ~ '^lsp plane. We set tan/3 = 2, sgn(/i) = — 

and = 100 TeV. It is greater than 0.7 and 0.9 in the red and blue regions, respectively. 


In Fig. 3d, we also vary the gluino mass as well as the LSP mass, showing the region where 
the branching fraction of the gluino into a Higgsino-like neutralino and a gluon is high. Here we 
set (/i/GeV,tan/3,m^^/TeV) = (—Mi — 100,2,100), with Mi being real and positive. We see. 
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rrig — m^o < 300 GeV and m^o — m^o ~ 100 GeV, can lead to efficient Z prodnction, BF{g —>■ 

gZxi) ^ 1 . 

3 LHC Signals 



Figure 4: Gonstraints at the 95% GL of the simplihed model. The black solid line shows the 
ATLAS multi-jets+MET constraints; the blue dashed line shows the GMS Z+MET constraint; and 
the green dotted line shows a combined ZZ constraint including the ATLAS four lepton, GMS four 
lepton, and GMS Z + dijet channels. The red regions show la and 2a parameter estimation from 
the ATLAS excess of the Z channel. 

The essential features of the present mass spectrum discussed in the last Section can be reduced 
into a simplihed model with the decay chain: g —)■ gx 2 and X 2 ^Xi- We assume m^o —m^o = 100 
GeV and take the branching fraction of each decay as 1 for simplicity. 

As discussed before, this simplihed model can produce sufficient amount of Z bosons, which we 
expect to be consistent with the ATLAS Z excess. On the other hand, the counterpart search by 
the GMS collaboration [2] should place severe constraints on it. In addition, this simplihed model 
can easily produce up to six jets due to the gluons produced by the gluino decay and the hadronic 
decay modes of Z bosons. Therefore, constraints by multi-jets+MET channels could potentially 
be important. Dilepton+Dijet+MET and four-lepton+MET analyses could also be relevant due to 
the (semi-)leptonic decay modes of Z bosons. Single lepton+MET channels are expected to be less 
important due to the second lepton veto. 

To study the htting region and various constraints, we generate this simplihed model with up to 
one extra parton in the matrix element using MADGRAPH 5 v2.1.2 [34] interfaced to Pythia 6.4 [35] 
and Delphes 3 [36] (which has Fast Jet incorporated [37]). The MLM matching [38] is applied with 


7 


















a scale parameter set to a quarter of the gluino mass. The parton distribution functions (PDFs) 
from CTEQ6L1 [39] are used. The gluino production cross sections are calculated at next-to-leading 
order (NLO) in the strong coupling constant, adding the resummation of soft gluon emission at 
next-to-leading-logarithmic accuracy (NLO+NLL) by using NLL-fast v2.1 [40]. 

In Fig. 4, we show Icr and 2cr parameter estimations from the ATLAS Z excess data [1]. For 
this fitting, we estimate the number of SUSY signal events for the ATLAS cut, which requires a 
same-flavor opposite-sign dilepton pair with its invariant mass in the Z mass range (81 GeV < 
mu < 101 GeV), two jets, and MET larger than 225 GeV. Additional cuts include the large scalar 
sum (i^T > 600 GeV) of the transverse momenta of all signal jets and the two leading leptons, and 
large azimuthal angular separations (A0(jeti,2, E™®®) > 0.4) between each of the leading two jets 
and the MET direction. With the observed number 16 (13) and the SM expectation value 4.2 ± 1.6 
(6.4 ± 2.2) for the dielectron (dimuon) channel, we construct a variable. The regions <2.3 
and 6.0, corresponding to 68th and 95th percentile of the Ghi-Squared distribution with two degrees 
of freedom, respectively, are referred to as Icr and 2cr fitting regions in the Figure. 

Similarly, we also estimate the signal strength of the GMS Z-I-MET data [2], ATLAS multi¬ 
jets (2-6 jets)-|-MET data [4] and four-lepton-|-MET data by ATLAS [9] and GMS [10] and GMS 
Z-|-dijet-|-MET [41]. Then we estimate the exclusion curves at the 95% confidence level, by using 
the CLg prescription. We use the SM background estimations and its uncertainties provided by each 
reference. Regarding the ATLAS multi-jets (2-6 jets)-|-MET searches, we choose a channel which is 
expected to give the most stringent constraints on each parameter point among 15 signal regions. We 
combine the four-lepton-|-MET data by ATLAS [9] and GMS [10] and GMS Z-|-dijet-|-MET [41] and 
show ZZ channel exclusion line in Fig. 4. Other constraints such as ATLAS large jet multiplicities 
(7-10 jets)-|-MET [15] and ATLAS Z-|-dijet-|-MET [42] are found less important. 

In Fig. 4, the relatively small mass splitting region [rrig — mNusp ^ 2m4) is of our true interest, 
because in this region the gluino decay branching fraction in the split MSSM can be very close to 1 
(see Fig. 3), and hence justifies our use of this simplified model. We see from the Figure that within 
the justified parameter region, there is a substantial parameter space that is consistent with ATLAS 
Z excess and not excluded by the various constraints, apart from the GMS Z-I-MET counterpart 
search. There is even a small parameter region consistent with both the ATLAS Z excess and the 
GMS counterpart exclusion limit. 

4 Summary and Discussion 

In this paper, we study the possibility of explaining the recent ATLAS Z excess in the MSSM 
spectrum. We study the gluino and neutralino decays in the case that sfermions are heavier than 
the gluino, assuming the Bino LSP and the Higgsino NLSP. We show that the small mass difference 
between the gluino and neutralinos and/or large stop mass can relatively enhance the radiative 
gluino decay g —)■ gH^. In this model, while the LSP is Bino-like neutralino, the direct gluino decay 
into the LSP is relatively suppressed. Motivated by this feature, we explored the simplified model 
to explain the ATLAS Z-I-MET excess. We found that the gluino mass around 800-1000 GeV and 
^NLSP ^ ~ 2mt can well explain the ATLAS Z excess without conflicting with other SUSY 

searches, apart from the GMS Z-I-MET search. In such a region the gluino radiative decays into X 2 3 
are dominated, as seen in Fig. 3, and the present simplified model well describes the realistic gluino 
decay chains. Therefore we can conclude the very simple MSSM spectrum may explain the ATLAS 


Z+MET excess. 

Moreover the present MSSM mass spectrum has another advantage. As pointed out in Ref. [12], 
the compressed mass spectrum can well fit the observed distributions of MET and the scalar 
sum of transverse momenta i/x in the signal region of the ATLAS Z+MET search. This feature is 
also the case for the present MSSM model. The mass spectrum of our interest, thus, can account 
for not only the number of the ATLAS Z excess without conflicting with the major SUSY searches 
but also the more detailed behaviours of the excess. 

However there is a subtlety when we consider the Z+MET search by the CMS, which seems to 
exclude a large portion of the best-£t region for the ATLAS Z+MET signals. Although the ATLAS 
and CMS searches are similar to each other, the ATLAS search makes relatively more account of the 
hadronic activity. This difference leads to slightly different LHC constraints between the ATLAS 
and CMS searches. Then there is a tiny region where the ATLAS excess can be explained and the 
CMS constraint is evaded. However the CMS exclusion limit and the ATLAS best fit region is very 
close, and it is hard to conclude the both searches can be consistent within this model, taking into 
account possible uncertainties of our fast simulations. A more detailed detector simulation will be 
needed to estimate more precise constraints. 

The LHC Run H will provide more obvious test for this model. The production cross section 
of the gluinos at 13 TeV LHC is enhanced by around 10, compared to the 8 TeV run, assuming 
the gluino mass is around 1 TeV. If the dilepton excess really comes from the SUSY particles, the 
number of the SUSY events for the integrated luminosity f Cdt = 10 fb“^ and the center-of-mass 
energy ^/s = 13 TeV will be around 100 with the same event selections as the present ATLAS 8 TeV 
dilepton plus MET search. The main SM background comes from tt production and its number is 
estimated to be around 50. Assuming the systematic uncertainty of the background estimation is 
30% as in the case of the 8 TeV result, a 5cr or more excess will be observed. Thus, the LHC Run 
H will provide a very clear test of the present SUSY model. 

It will be worth noting that this mass spectrum may provide a good Bino-like dark matter can¬ 
didate. The abundance and detection of the dark matter is quite sensitive to the other parameters, 
such as the Wino mass, CP phases and tan/?, and its detailed study is beyond the scope of the 
present paper. 

In this study, we assume the split SUSY-like spectrum in which the scalar tops play a dominant 
role in the gluino decay, and the Bino and Higgsino are the LSP and NLSP, respectively. Although 
this assumption is well motivated, it is interesting to investigate more generic types of MSSM spectra, 
and it will be done elsewhere. 
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